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We tackle the problem of unraveling the algebraic structure of computations of effective Hamil- 
. tonians. This is an important subject in view of applications to chemistry, solid state physics or 

quantum field theory. We show, among other things, that the correct framework for these compu- 
tations is provided by the hyperoctahedral group algebras. We define several structures on these 
algebras and give various applications. For example, we show that the adiabatic evolution opera- 
tor (in the time-dependent interaction representation of an effective Hamiltonian) can be written 
, naturally as a Picard-type series and has a natural exponential expansion. 

> 



Abstract 



^ ■ Introduction 



We start with a short overview of the classical theory of Chen calculus, that is, iterated integral com- 
putations. The subject is classical but is rarely presented from the suitable theoretical perspective -that 
is, emphasizing the role of the convolution product on the direct sum of the symmetric groups group 
', algebras. We give therefore a brief account of the theory that takes into account this point of view 

-this will be useful later in the article. Then, we recall the construction of effective Hamiltonians in the 
time-dependent interaction representation. 
' The third section is devoted to the investigation of the structure of the hyperoctahedral group algebras. 

Although we are really interested in the applications of these objects to the study of effective Hamiltonians, 
and although the definitions we introduce are motivated by the behavior of the iterated integrals showing 
up in this setting, we postpone the description of the way the two theories interact to a later stage of 
the article. Roughly stated, we show that the descent algebra approach to Lie calculus, as emphasized in 
Reutenauer's can be lifted to the hyperoctahedral setting. This extends previous works [T^ITlfSllSlfTni 
on the subject and shows that these results (focussing largely on Solomon's algebras of hyperoctahedral 
groups and other wreath product group algebras) are naturally connected to the study of physical systems 
through the properties of their Hamiltonians and of the corresponding differential equations, very much 
as the classical theory of free Lie algebras relates naturally to the study of differential equations and 
topological groups. Notice however that the statistics we introduce here on hyperoctahedral groups seems 
to be new -and is different from the statistics naturally associated to the noncommutative representation 
theoretic approach to hyperoctahedral groups, as it appears in these works. 

The fourth section studies the effective adiabatic evolution operator and shows that it can be expanded 
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as a generalized Picard series by means of the statistics introduced on hyperoctahedral group^. As a 
corollary, we derive in the last section an exponential expansion for the evolution operator. Such expan- 
sions are particularly useful in view of numerical computations, since they usually lead to approximating 
series converging much faster than the ones obtained from the Picard series. 

Acknowledgements. We thank Kurusch Ebrahimi-Fard, Jean- Yves Thibon and an anonymous referee. 
Their comments and suggestions resulted into several improvements with respect to a preliminary version 
of this article. 



1 The algebra of iterated integrals 

Let us recall the basis of Chen's iterated integrals calculus, starting with a first order linear differential 
equation (with, say, operator or matrix coefficients): 

A'{t) = H{t)A{t), A{0) = 1 

The solution can be expanded as the Picard series: 

t t ti 

A{t) = l + j H{x)dx + j j H{ti)H{t2)dtidt2 + ... + J H{ti)...H{tn) + ... 

AJ, 

where := {0 < t„ < ... < ti < t} and where the measure dti . . .dtn is implicit. Notice, for further 
use, that when we allow for a general initial condition A{x) = 1, with possibly x = — oo, the integration 
simplex should be replaced by A[f'*' :— {x < t„ < ... < ti < t\. 

Solving for A{t) = exp(r2(i)) (see O [H]), and more generally any computation with A{t), requires 
the computation of products of iterated integrals of the form: 



< cr >:= j if(t^(i))...i/(t^(„)), a = (ct(1), cr(n)) e 5„, 
A* 

where Sn stands for the symmetric group of order n. Notice that we represent an element a in S'„ by the 
sequence {cf{\), ...,a{n)). 

In general, for any H — J2 S fJ-a ■ ^ ^ '.^ ®Q[Sn], the direct sum of the group algebras of the 

n a^Sn n 

symmetric groups Sn over the rationals, we will write < fi > for ^ ^ fig-- < a >. This allows, for 
example, to write A{t) as < / >, where / := "^{1, n) is the formal sum of the identity elements in the 

n 

symmetric group algebras. 

The formula for the product of < cr > with < /? > is a variant of Chen's formula for the product of 
two iterated integrals of functions or of differential forms (a proof of the formula will be given in Section [3] 
in a more general framework; the formula also holds for arbitrary integration bounds, that is when A^ is 
replaced by aI^ '*'): 

< a > ■ < f3 >=< a * (3 >, 

where a * f3 is the convolution produdU of the two permutations, that is, for a E Sn, P & Sm'. cr * f3 
is the sum of the ("^'") permutations 7 G Sn+m with st{j{l), ...,7(71)) = (a(l), ...,cr(ri)) and st{j{n + 
1), ...,7(n -|- m)) = (/3(1), ...,P{m)). Here, st stands for the standardization map, the action of which on 
sequences is obtained by replacing (ii, i„), e N* by the (necessarily unique) permutation a G Sn, 
such that cr{p) < a{q) for p < q\i and only if ip < iq. In words, each number ij is replaced by the position 



^Picard series are often referred to as Dyson or Dyson-Chen series in the literature, especially in contemporary physics, 
but we prefer to stick to the most classical terminology 

■^This is one possible definition of the convolution product, there are several equivalent ones that can be obtained 
using the various natural set automorphisms of the symmetric groups (such as inversion or conjugacy by the element of 
maximal length). They result into various (but essentially equivalent) associative algebra structures on the direct sum of 
the symmetric groups group algebras, see e.g. jl3) 
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of ij in the increasing ordering oi ii, . . . , in- If we take the example of (5, 8, 2), the position of 5, 8 and 
2 in the ordering 2 < 5 < 8 is 2, 3 and 1. Thus, st(5, 8, 2) = (2, 3, 1). For instance, 

(2,3,1)*(1) = (2, 3, 1,4) + (2, 4, 1,3) + (3, 4, 1,2) + (3, 4, 2,1), 

(1,2)* (2,1) = (1,2,4, 3) + (1,3,4,2) + (1,4, 3,2) + (2,3,4,1) + (2, 4, 3,1) + (3,4,2,1). 

The convolution product relates to the shuffle product 111 [H] as it appears in Chen's work T and 
in the parametrization of the product of iterated integrals of functions or differential forms. Indeed, for 
cr, P as above: 

a-Up-^ ^ (am/3[n])-i 

where /3[ti] stands for the sequence {(3{l) + n, /3{m) +n). Recall, for completeness sake, that the shuffle 
product A[1]B of two words (or sequences) A = aA' ~ aa2...ak, B = hB' = 662. ..6; is defined recursively 

by 

AWIB = a(yl'LUB) + h{A[hB') 

we refer to [U] and [IS] for details. Associativity of * follows immediately from the definition, the unit 
is 1 e So = Q, and the graduation on S = 0Q[S'n] is compatible with *, so that: 

n 

Lemma 1.1. The convolution product provides S with the structure of a graded connected associative 
(hut noncommutative) unital algebra. 

For completeness, recall that connected means simply that Sq — Q. In fact, S carries the richer 
structure of a Hopf algebra, and as such is referred to as the Malvenuto-Reutenauer Hopf algebra [T^. 
From the point of view of the theory of noncommutative symmetric functions, the elements of S should 
be understood as free quasisymmetric functions [8J. This definition of the convolution product on S al- 
lows, for example, to express simply the coefficients of the continuous Baker-Campbell-Hausdorff formula 
(compare with the original solution '16'): 

m =< log(^) > • 

Here log(/) identifies, in S, with the formal sum of Solomon's (also called Eulerian or canonical) idem- 
potents pS] . We refer to [5D1 dSl [HJ [551 [S] for an explanation and a Hopf algebraic approach to these 
idempotents and, more generally, for a Hopf algebraic approach to Lie computations. We will return 
later with more details to Solomon's idempotent but mention only, for the time being, that one of the 
main purposes of the present article is to extend these ideas to the more general framework required by 
the study of effective Hamiltonians. 

2 Iterated integrals in time-dependent perturbation theory 

The ultimate aim of quantum physics is the knowledge of the eigenvalues and eigenstates of the Hamil- 
tonian H describing a physical system. In most cases, the eigenvalue problem cannot be solved exactly, 
but the eigenvalues and eigenstates of a simpler and closely related Hamiltonian Hq are known, at least 
numerically. Then, H can be rewritten H = Hq + V , where V is referred to as the perturbation term. For 
example, in molecular physics, Hq is the Hamiltonian describing the interactions of A^e electrons with Nn 
nuclei, whereas H describes the interaction of the electrons with themselves and with the iV„ nuclei. 
In that case, the perturbation V describes the electron-electron interaction, this simple approach paving 
the way to most of the numerical methods in the field. 

In other terms, perturbation theory provides a systematic way to calculate an eigenstate of H from 
an eigenstate of Hq. In the time-dependent perturbation theory, we first define a time-dependent Hamil- 
tonian H{t) = Hq + e~^l*ly. When e is small, this means physically that the interaction is very slowly 
switched on from t = —00 where H{—oo) = Hq to i = where H{Q) = H. The basic idea is that, if 
e is small enough, then an eigenstate of Hq can be transformed into an eigenstate of H by the time- 
dependent perturbation e~'^^*^V. For example, the ground state |<i>o > (the eigenstate associated to the 
lowest eigenvalue Eq of Hq) should hopefully be transformed into the ground state of H. We also assume 
for the time being that the ground state is non degenerate, that is that the eigenspace associated to the 
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highest eigenvalue Eq is one-dimensional. We write from now on |$i > for the eigenvectors of Hq and 
assume that the eigenvalues Ei are ordered increasingly {Ei < E'i+i). 

To implement this picture of perturbation theory, the time-dependent Schrodinger equation id\'^ s{t)) / dt ~ 
H{t)\'i> s{t)) should be solved. However, the solutions |5'5(t)) of this equation vary like e~*^''*|$o} for 
large negative t (where |$o} satisfies iJol^'o) = -^'ol'I'o)) and have no limit for i — » — oo. To compensate for 
this time variation, one looks instead for = e'^°*|*s(i)) that satisfies id\^{t))/dt = Hi^{t)\^ {t)) , 

with iJie — e'^°*V^e~*^''*e~^l*l. Now Hi^{—oo) = 0, and oo)) makes sense. Using Hi^, we can 
start consistently from the ground state |$o) of Hq and solve the time-dependent Schrodinger equation 
with the boundary condition cx))) = |$o)- When no eigenvalue crossing takes place, |$o} should be 
transformed into the ground state |5'(0)} of H . 

At first sight, solving the time-dependent Schrodinger equation for \^{t)) does not look simpler than 
solving the time-independent Schrodinger equation with the Hamiltonian H . However, if V is small 
enough, fairly accurate approximations of the true eigenstates can be obtained from the first terms of 
the perturbative expansion of {^{t)). In general, instead of calculating directly \^{t)), it is convenient to 
define the unitary operator Uc{t) as the solution of idUe{t)/dt = Hi^{t)Ue{t), with the boundary condition 
Ue{—oo) = 1. Thus, |5'(t)} = Uf{t)\^o). The connection with iterated integrals appears when solving 
iteratively the equation for Ue{t): 

US) = ! + (-*)/ AhHu{h) + {-if I dti r AhHi,{h)Hu{t2) + . . . 

J — OO J —OO J —oc 



A straightforward calculation [TT] gives us 



um^) It M V V |$.j(j>.,jT/|$.„_,)...(^jy|^.,)($jF|$o) 

n— 1 . . .?„ 

where we use the completeness relation 1 = |<i>i}(<i>,i|. However, it immediately appears that this 
expression has no limit for e — > because the sum over all ip contains the terms ip ~ 0, for which the 
denominator has a factor pie. In 1951, Gell-Mann and Low [TU] conjectured that 

*gl) = hni' 



.0 ($o|f/.(0)|$o) 



exists and is an eigenstate of H . The fact that the Gell-Mann and Low formula indeed removes all the 
divergences of equation ([T]) was proved much later by Nenciu and Rasche [IH] . 

The above scheme works nicely when the ground state of Hq is non degenerate. When it is degenerate, 
that is when the eigenspace Eq associated to the lowest eigenvalue of Hq has dimension > 1, the problem 
is more subtle, see |17l O I15| . To understand the algebraic phenomena underlying time-dependent 
perturbation theory for degenerate systems is actually the purpose of the present article. 

Let us write P for the projection on this eigenspace. The natural extension of the Gell-Mann and 
Low formula then reads as a definition of a "Gell-Mann and Low" operator acting on the degenerate 
eigenspace Eq: 

Ugl ■■= lim?7e(0)P(Pf/,(0)P)-\ (2) 

£— >0 



or, Ugl = limf7GL(e), UgU^) U,{0)P{PUe{0)Py^ . It can be shown that the operator PU^[Q)P is 

invertible within the image of P if no vector in the image of Ue{0)P is annihilated by P [1]. We will 
assume that this property is satisfied by the systems we consider. The operator Ugl was first proposed 
by Morita in 1963 [I?]. It shows up e.g. in the time-dependent interaction representation of the effective 
Hamiltonian Hcs := lim PiJ/7e(0)P[Pt/(;(0)P]~^ classically used to solve the eigenvalue problem. 

This operator Ugl is the one we will be interested in here, postponing to further work the analysis of 
concrete applications to the study of degenerate systems. In other terms, we will investigate and unravel 
the fine algebraic structure of the iterative expansions of Uf{t) and UGL{i)- 
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3 Wreath product convolution algebras 



Let us explain further our motivation. In the previous section, we observed that the study of effective 
Hamiltonians leads to the study of Picard-type expansions involving the operators Hi^{t) and PHi^{t) or, 
equivalently, A{t) — i(l — P)Hi^{t) and B{t) := iPHi^it). Expanding these expressions will lead to the 
study of iterated integrals involving the two operators A{t) and B{t) such as, say: / A{t2)B{t-i)A{ti). 

A* 

The idea underlying the forthcoming algebraic constructions is to encode such an expression by a signed 
permutation and to lift computations with iterated integrals to an abstract algebraic setting: in the 
previous example, the signed permutation would be (2, 3, 1) (see below for precise definitions). 

In more abstract (but equivalent) terms, iterated integrals on two operators are conveniently encoded 
by elements of the hyperoctahedral groups, whereas iterated integrals on k operators would be conve- 
niently encoded by more general colored permutations or elements of the wreath product of 5„ with the 
cyclic group of order k. Recall the definition of the hyperoctahedral group i?„ of order n. The hyper- 
octahedral group is the group defined either as the wreath product of the symmetric group of order n 
with the cyclic group of order 2, or, in a more concrete way, as the group of "signed permutations" the 
elements of which are written as sequences of integers i G N* and of integers with an upper bar i, i G N*, 
so that, when the bars are erased, one recovers the expression of a permutation. The composition rule 
is the usual one for permutations, together with the sign rule for bars: for example, if a G i^a = (2, 3, 1) 
and /3 = (3,1,2), then: 

/3oct(2) = /3(3) = 2 2, 
^oa(3)= /?(!) = 3. 

By analogy with S, we equip B := ^i?„ with the structure of a graded connected (associative but 

n 

noncommutative) algebra with a unit. This algebra structure agrees with the ones introduced in [Tl ll9l [^ 
(possibly up to an isomorphism: for example, the relationship between the product we consider and 
the shifted shuffle product of [19] reflects the relationship between the convolution and (shifted) shuffle 
product on S that we recalled in the first section of the present article) . 

The standardization st of a signed sequence w (i.e. a sequence of integers and of integers marked 
with an upper bar) is defined analogously to the classical standardization, except for the fact that upper 
bars are left unchanged (or, equivalently, have to be reintroduced at their initial positions after the 
standardization of the sequence w has been performed, where we write w for w where the upper bars 
have been erased). For example, si(2, 7, T, 2) = (2, 4, T, 3). 

Definition 3.1. Let a, (3 belong to i?„, resp. Bm. Their convolution product is defined by: 

a * /3 :— f 

T 

where f runs over the ("^'") elements of Bn+m with si(f (1), f{n)) = a, st{f{n + 1), f{n + m)) = /3. 
For instance, 

(2,3,l)*(i) = (2, 3, 1, 4) + (2, 4, 1, 3) + (3, 4,1, 2) -f (3, 4, 2,1), 

(1, 2) * (2, 1) = (1, 2, 4, 3) + (1, 3, 4, 2) + (1, 4, 3, 2) + (2, 3, 4, 1) + (2, 4, 3, 1) + (3, 4, 2, 1). 

Notice that this definition is dictated, for us, by iterated integrals computations, similarly to the classical 
one-Hamiltonian case dealt with in the first section. Indeed, let A{t),B{t) be two time-dependent oper- 
ators. For a G Bn, let us write < a > for the iterated integrals obtained by the usual process, with the 
extra prescription that upper indices (empty set or bar) in a indicate that the operator used at the cor- 
responding level of the integral is A or B, so that e.g., a = (3, 1, 2) is associated to: / B{t3)A{ti)B{t2). 

For an arbitrary 7 = S as • f G B, we write < 7 > for ^ ^ Og- < a > . 

Proposition 3.2. The product of two iterated integrals < a > x < f3 > is given by: 

< a > X < (3 >=< CT * /3 > 
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This formula is a noncommutative variant of the classical Chen formulas for the product of iterated 
integrals of differential forms [7^ . It includes as a particular case the formula for the product of two iterated 
integrals depending on a single time-dependent Hamiltonian given in the first section of the article. We 
detail the proof for the sake of completeness, and since the formula is crucial for our purposes. 

For a permutation a we denote by a the same permutation without bars (e.g. if ct = (2, 3, 1), then 
(J = (2,3,1)) and we define X(tcr(i)) = ^(^cr(i)) if 5-(i) has no bar and X{t„i^i'j) = i?(tcr(i)) if ^(i) has a 
bar. Therefore, 



< a > X < f3 > = / dti . . . dtnX{t„(^i)) . . . X{t„(n)) / dtn+l . ■ ■ dtn+jnX{tn+p{l)) ■ ■ ■ X{tn+p(,n))- 

Ja Jo Jo Ja 

By Fubini's theorem, this can be rewritten as the integral of X{t^(^i-^) . . .^(i„+/3(m)) over the domain 
X A^j. The idea is now to rewrite this domain as a sum of ("^j"') domains isomorphic to A*j^„j. For 
instance, the product of the domain < tn < ■ ■ ■ < ti < t with the domain < t„+i < t is the sum of the 
n-\-l domains obtained by inserting t„+i between and ^i, then between ti and t2, up to between i„ and 
t. More generally the product of A^ by is the sum of all the domains obtained by "mixing" the two 
conditions < tn < ■ ■ ■ < ti < t and < tn+m < • ■ ■ < ^n+i < t, i.e. by ordering the n + m variables U 
so that these conditions are satisfied. If p{i) is the position of variable ti in one of these orderings (where 
the variables are ordered from the largest to the smallest), the conditions imply that p{l) < ■ • • < p{n) 
and p{n + 1) < • • • < p{n + m). For example, if < ^2 < < ^ and < t4 < t^ < t, for the domain 
^ t^ < t2 < ti < t^ < t, ts is in the first place (i.e. largest), ti in the second, t2 in the third and t4 in 
the fourth (smallest), and the permutation is p = (2, 3, 1, 4). In general, we get: 

A^ X A^ = U{(^-(1)' ^r(«+m))|0 < t„+™ < ... <ti< t}, 

T 

where r runs over the permutations in Sn+m such that t(1) < • • ■ < r(n) and T{n + 1) < ■ • ■ < 
T(n + m). Equivalently, r runs over the permutations such that: st(r(l), . . . ,T(n)) = (1, . . . ,n) and 
st{T{n + 1), . . . , T(n + m)) = (1, . . . , to). Now, 

/ X{x„(^i-))...X{Xcr{n))X{Xn+(3{l)) ■ ■ ■ X {Xn+f3(^rn)) 

{{xi=t^fiy...,Xn + m = t^^ri + m))\0<t„ + m<---<tl<t} 
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X{tT{a{l)))---X{tT{a{n)))X{tT-(^n+f3{l)))---X{tT(^n+l3{m))) 

{0<t„+m<...<ti<t} 

SO that finally, taking into account the bars of the permutations (that is the fact that X is A or B, 
depending only on its position in the sequence X{tr(a(iy))-X{tr{a{7i)})X{tr(n+i3(iy))-X{tr(n+i3(m)))), we 
obtain < a > x < f3 >= J2j < 7 >j with 5^(7(1) . . . 7(«)) = (5'(1), ■ ■ • , ^(n)) and st{^{n + 1) . . . ^{n + 
to)) = (/3(1), . . . , /3(n)). This concludes the proof. 

Remark 3.3. The same proof leads to a general noncommutative Chen formula. Let Ai, . . . ,An and 
Si, ... , Bm be noncommutative (e.g. matrix-valued) functions and let 

Aa = dti . . . dtnAi{ta{l)) ■ ■ . An{ta(n)), 

Ja Jo 

Bfj = I dti... I dtm Bi{t 13(^1)) ... Brn{t[3{rn)), 

Jo Jo 

pt rt„ + m-l 

{AB)a = / dti... dt„+„iAi(to.(i)) . . . y4„(to.(„))i?i(ic,(„_,_i)) . . . i3,„(tcr(,i+m)), 

Jo Jo 

then Aa x Bp = {AB)a*fj. 

Proposition 3.4. The convolution product provides B with the structure of an associative (but noncom- 
mutative) algebra with a unit. 
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The Proposition can be checked directly from the combinatorial definition of the convolution product, 
we refer to the original proofs [1, 19J. In our setting, it also follows from the associativity of the product 
of iterated integrals. Notice that the general noncommutative Chen formula would relate similarly to the 
algebraic structures on colored permutations and wreath product group algebras introduced in [19j . 

4 Progressions and regressions 

Modern noncommutative representation theory originates largely in the work of Solomon [29J . From this 
point of view, it is natural to partition hyperoctahedral groups into "descent classes" , similarly to the 
partition of symmetric groups into descent classes (such a partition is often referred to as a statistics on 

Sn). 

Recall that a permutation a G Sn has a descent in position i < n if and only if a{i) > a{i + 1). The 
descent set Desc{a) of a is the set of aX\ i < n such that a has a descent in position i. The partition into 
descent classes read: Sn ~ U Wi Desc{a) — /}. The descent algebra T) is the linear span of Solomon's 

/C[ri-1] 

elements D'g :— ^ cr, where S" C [n — 1] and n e N* (with the convention = 1). It is 

(j^S„,Desc{a)(ZS 

provided with a free associative algebra structure by the convolution product * on S D P, see [351 Chap. 9]. 
This algebra has various natural generating families as a free associative algebra -for instance, the family 
of the . It is therefore also isomorphic to the algebra of noncommutative symmetric functions Sym, 
from which it follows that the structure theorems for these functions can be carried back to the descent 
algebra -a point of view introduced and developed in ^ and a subsequent series of articles starting with 

m- ' ' 

The corresponding descent statistics on Bn is obtained by considering the total order n < n — \ < 
... < 1 < 1 < ... < n. A signed permutation a G Bn has a descent in position i < n if and only 
if a{i) > a{i + 1) [Ml Def. 3.2]. Descent classes are defined accordingly. The problem with this 
noncommutative representation theoretical statistics and with the corresponding algebraic structures is 
that they do not fit the needs of iterated integral computations for effective Hamiltonians, as we shall 
see in the forthcoming sections. Neither do the generalized descent algebras of [T^ Sect 5]. Notice 
that this is not the case when symmetric groups are considered: the statistics of descent classes fits the 
needs of noncommutative representation theory as well as the needs of Lie theoretical computations, as 
emphasized in [5S1 [S] . 

For this reason, we introduce another statistics on Bn. It seems to be new, and has surprisingly nice 
properties, in that it allows to generalize very naturally many algebraic properties of symmetric groups 
descent classes. 

We say that an element a — (q;(1), a{n)) G Bn has a progression in position i if either: 

1. \a{i)\ < \a{i + 1)| and a{i + 1) G N* 

2. \a{i)\ > \a{i + 1)| and a{i + 1) G N* 

Else, we say that a has a regression in position i. Here, the operation | | is the operation of forgetting 
the bars, so that e.g. |6| = 6. The terminology is motivated by the quantum physical idea that particles 
(associated to unmarked integers) propagate forward in time, whereas holes (associated to marked integers 
in our framework) propagate backward. We refer the reader to Goldstone diagrams expansions [TT] of 
the Gell-Mann Low eigenstate I^'gl > for further insights into the physical motivations. Further details 
on these topics are contained in the following sections of this article, but we do not develop here fully the 
physical implications of our approach, the focus being on their mathematical background. 

We write Reg{a) for the set of regressions of a. For example: Reg{A, 3, 5, 6, 2, 1) = {2, 5} since the 
sequence (4, 3, 5, 6, 2, 1) has only two regressions, in positions 2 and 5. For an arbitrary subset 5 of [n — 1], 
we mimic now the descent statistics and write Kg :— ^ a. It is also convenient to introduce 

a£B„,B£g{<7)=S 

the elements Tg -.^ J2 ^ = Ru- 

aeB„,Reg{a)CS UCS 

Lemma 4.1. The elements Rg (resp. Tg), S (^[n — 1], form a family of linearly independent elements 
in the group algebra Q[i?„]. 
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The first assertion follows from the very definition of the Rg, since it is easily checked that {a G 
Bn, Reg{a) = S*} ^ for any 5 C [n — 1]. The second case follows from the Mobius inversion formula: 

u<zs 

where \S\ stands for the number of elements in S. 
Lemma 4.2. We have, for 5 C [n — 1], [/ C [m — 1]; 

Tin , rjim rji 
S*-^U — -'-SU{n}U{U+n), 

where U + n — {u + n,u €z U}. 

Indeed, by definition, for a e -B„, /3 G -B™, with Reg{d) — X <Z S, Reg{(3) — Y C U , a * (3 = 

T _ 

where f runs over the elements of Bn+m with si(f(l), ...,f(n)) = tr and st{f{n + 1), ...,f(n + m)) ~ (3. 
In particular, for any such f and by definition of the standardization process: 

Reg{f) CXU{n}U(y + n). 

Conversely, any f G Bn+m appears in the expansion of st(f(l), T(n)) * st{f(n + 1), f(n + m)) 
by the very definition of * and does not appear in the expansion of any other product a * (3 with 
Reg(a) = Reg{st{f(l), f (n))), Reg{0) ~ Reg{st{f{n+ 1), f{n + m)), from which the lemma follows. 

Corollary 4.3. For S, U as above, with the notation of the previous sections: 



so that: 



Theorem 4.1. T/ie linear span TZ of the elements Tg (equivalently, of the R^), n G N, S" C [n — 1], is 
closed under the convolution product in B. This algebra, referred to from now on as the (hyperoctahe- 
dral) Regression algebra, is isomorphic to the descent algebra T> and to the algebra of noncommutative 
symmetric functions Sym. 

The second part of the Theorem follows from the product rule in T), that reads: 

Ug* Ujj - i^5u{n}u(;7+n)- 

The proof for this last identity can be obtained similarly to the one in Lemma W?^ -see also [25] . 

Now we study in more detail the elements R^ that will play an important role in the following. The 
lowest order i?^ are 

Rl = (l,2) + (i,2) + (2,i) + (2,l), 

Rl = (1, 2, 3) + (1, 2, 3) + (1, 3, 2) + (1, 3, 2) + (2, 1, 3) + (2, 1, 3) + (2, 3, 1) + (2, 3, 1) 
+ (3, 1,2) + (3, 1, 2) + (3, 2,1) + (3, 2,1). 

We first observe that, if (i G Bn is a term of R^ (and therefore has no regression), then the barred 
integers of a are entirely determined by the permutation a = (|a'(l)|, . . . , |(7(n)|), except for <t(1). Indeed, 
by definition of a progression, a{i + 1) G N* if a{i) < a{i + 1) and a(i + 1) G N* if cr(i) > a(i + 1). In 
other words, a{i + 1) G N* iff cr has a descent at i. The integer (7(1) is not determined by a and can be 
barred or not. Therefore, the number of terms of is 2 • n!. 
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5 A Picard-type hyperoctahedral expansion 



When it comes to expand ^gl or Ugl, as introduced in Section 2, the classical strategy introduced by 
Goldstone (at least for nondegenerate states, that is for "^gl [H]) consists in appealing to the hole/particle 
duality of quantum physics. Goldstone's theory was generalized to degenerate states by Michels and 
Suttorp [15] , but this part of the theory has remained largely in infancy and relies on shaky mathematical 
grounds. The purpose of this section is to show that hyperoctahedral groups provide a convenient way 
to derive and study such expansions, so as to build the foundations of a group-theoretic approach to the 
perturbative computation of the ground states of physical systems, with a particular view toward the 
degenerate case. 

To sum up, we want to compute Ugl{£) — U^{0)P{PU^{0)P)^^ . Let us write Hi for —iHi^ and 
A{t) (1 - P)H,{t), B{t) := -PH^{t) (notice the -1 sign in the definition of B). From the Picard 
expansion, we have: 

tl 

UM = ^+ j H,{x)dx+ j J H,{h)H,{t2)dtidt2 + ...+ J H,{h)...H,{tn) + - 

— oo — oo — oo ^[-oo>0] 

In this section (following a suggestion by the referee whose remarks helped us to simplify notably the 
presentation of the following computations -we take the opportunity to thank him or her warmly) , we 
introduce a new encoding of iterated integrals involving A, B and Hi. 
The notation is best explained through an example: 

tl t2 ts ^4 

J A{t3)H,{ti)B{t2)B{t5)H,{ti)dhdt2dt3dUdt5 =: [2,3,1,5,4]. 

— oo — oo — oo — oo — oo 

Concretely, in an arbitrary iterated integral involving A, B and Hi, we look recursively at the positions 
ii, ...,in of ti,...,tn in the integrand and decorate ij with a bar (resp. a hat, resp. no decoration) if the 
corresponding operator is B (resp. Hi, resp. A). In our example, ti (resp. t2---) is in position 2 (resp. 
3...) in the product A{t-3)Hi{ti)B{t2)B{t^)Hi{t4) and appears as a parameter for Hi (resp. B), so that we 
map 1 to 2, 2 to 3, and so on. The so-obtained sequence of decorated integers is written inside brackets to 
avoid confusion with our previous notation. Notice that, if G -B„, < d >= [a~^{l), ...,a~^{n)] (where 
we use the definition of Section [3] for < ct >). We will also use some self-explaining multilinear extensions 
of this notation, so that, for example: [3, (1, 2) + (2, 1)] = [3, 1, 2] + [3, 2, 1], and so on. Notice in particular 
that the identity Hi ^ A — B translates formally into k = k — k so that, for example, [1, 2] = [1, 2] — [1, 2]. 

Theorem 5.1. The effective adiabatic evolution operator Ugl has the hyperoctahedral Picard-type ex- 
pansion: 

f/GL = limF+(l-P)( ^ <iq>)P 

rieN* 

Indeed, let us expand [1, ...,n] = J Hi{ti)...Hi{tn) with the A and B operators. In order to do 

^[-00,01 

so, we introduce still another notation: for a G Bk, k < n, we set: 



<a;n-k>= / X{t^(^i))...X{t^(^k))H,{tk+i)...H,{tn) 



where a[, x A^J^'*"**"' is a shortcut for: 

{{tl, ...,tn)\ - oo <tk < ... <ti <0, -oo < t„ < ... < tk+l < tcr(fc)}; 

where a stands, as usual, for the image of a in Sk (obtained by forgetting the decorations), and where 
= ^(^(T(i)) if "■(*) = and B{t^(^i^) else. For example, 

<(213);2>= J A{t2)B{ti)A{t3)H,{ti)H,it5), 
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<(231);2>= J A{t2)Bm)A{h)H,{U)H,{t^). 

-CX3<t3<t2<tl<0, -00<t5<t4<ti 

The more general symbols < X]n — k > are defined, as usual, by extending linearly these conventions to 
arbitrary elements X G k < n. 

Lemma 5.1. The symbol < a;n — k > can be expanded as: 

<a;n-k>^ [a'^l), ...,a-\i), {a-\i + 1), ...,a-\k))[H{k + l, ...,n)], 
where a G Bk and i :— a(k). 
Indeed: 

{(ti,...,i„)| -oo < ife < ... <<i <0, -oo<t„ < ...<tk+i <ta(fe)} = 

{(tl, ...,t„)| - oo < ifc < ... < io.(fe), -oo <tn< ... < tk+1 < tcy(k),ta{k) < ■■■ < <1 < 0} 
= {(il,i2, ■■■,tn)\ - OO < i^(fe) < ... < il < 0, (ta(/c) + l, ■■■,tn) G 
...,Ufe-<T(fe))| - OO < U„(k) < ... < Ul < ta(fc)} X {(Wi, ...,'y„_fe)| - oo < Vn-k < ■■■ < Vi < t„(^k)}}, 

where x stands for the cartesian product. Since the Cartesian product of simplices is reflected in the 
shuffle product (see e.g. the proof of Prop. . the Lemma follows. 
We then have: 

[i, ...n] = [1, 2, ...n] - [1, 2, n] = -[1, 2, n] + [1, 2, 3, ...n] - [1, 2, 3, n] 
= P[i,...h] + [1,2,3, ...n] - [1,2,3, 

= + [1,2,3,. ..n]- [1,2,3,. ..,n] + (-[2,(l)m(3,...,n)] + [2,(l)m(3,...,n)]). 

Now, from the recursive definition of the shuffle product: 

[1,2,3, ...,«] + [2, (l)m(3, n)] = [(l)m(2, 3,..., n)] 



A{t)dt] J B{h)H,{t2)...H,{tn-l) 

= -(l-P)<i?]>P[i,...,rr^l]. 
Here, we have used that P is a projection, so that (1 — P)B{t) = —(1 — P)PHi{t) — 0, to rewrite 



[ j A{t)dt] = {l-P)[j {A{t) + B{t))dt] = {l-P)<Rl>. 

— oo — oo 

We get finally, since < (l,2);n— 2>= [1, 2, 3, n] and (according to the Lemma [??T|) < (2, 1); n— 2 > = 
[2,(l)m(3,...,n)]: 

[i,...,n] =P[i,...,n] + (l-P) <i?J >P[i,...,rr^l] + (l-P) <i?^;n-2>, 

where the last identity follows, once again, from (1 — P)B{t) = (we won't comment any more on this 
rewriting trick from now on). 

The proof of the Theorem can be obtained along these principles by recursion. Let us indeed assume 
for a while that: 

< R^;n- k >=< > P[i, ...,n^]+ < R^'^^^n- fc - 1 > . 
Then we get, by induction: 

[i,...,n] =P[i,...,n] + (l-P) <Pj>P[i,...,n^l] + 
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(1 - F) < Rl > p[i, ...,n - 2] + ... + (1 - p) < > P[i] + (1 - P) < i?^ > . 

Since /7e(0) = X^nflj this implies 

oo 

UM = pc/e(0) + (i-p)^ <p^' > (p(c/,(o)-i) + i), 

or, since P^ = P: 

CX) 

n=l 

and the Theorem follows. 

So, let us check that the formula for < R^\n — k > holds. Let us consider an arbitrary element a £ Bk 
with Reg{a) — 0. Then, with the notation of Lemma l5. II 

<a;n-k>^ [^-^(l), ^-^(i), (^-^(i + 1), ^-^(fc)) 111 (fc + 1, fcT2, n)] 

-[5-1(1), (a-i(z + 1), ...,a-i(fc))m(fcTT, ffi, n)] 

Let us denote the first term by Ti and the second by T2, so that < a\n — k >= Ti — T2. To calculate Ti, 
let us use another (equivalent, the equivalence follows from the recursive definition of the shuffle product 
and is left to the reader) recursive definition of the shuffle product, namely: 

k 

aia2---ak\l^hh2...hi = ^ ai...ai6(ai+i...afeLLI52...6;). 

k 

Ti - + 1, {a-\j + 1), ...,a-i(fc))m(fcT2, n)] 

k 

= < /3j;n - fc - 1 > 

(a-"i(l), ...,5--i(j),A: + l,o--i(j + 1), ...,o--i(fc))-^ We notice then that, for / < k, 

(3,(1) < I3j{l + 1) ^ a{l) < ail + 1). 

In particular, /3j has no regression in position less than k. Now, /3j(fc + 1) = /3j(fc + 1) = J > I3j{k) = i, 
which implies that Pj has no regression in position k. Finally, f3j has no regression. 

Let us enumerate the number of (necessarily distinct) signed permutations f3j obtained in that way. 
There are 2(fc — 1)! elements a of i?^ with a given value j of cr(fc), where j runs from 1 to k. For a given 
a with a{k) — j, Ti provides k — j + 1 elements of R0~^^- Thus, the expansion of Ti provides (fc + 1)! 
different elements of R'^'^^ when a runs over Pg. 

The term T2 can be computed similarly. Using the same recursive formula for the shuffle product as 
above, we get: 

k 

r2 = 5][a-i(l),...,a-i(j),Fn, (a-i(j + l),...,a-i(fc))m(fcT2,...,n)] 

= [{<j-\l),...,<j~\k))[n(kTl,...,n)] 
- ^[(5-1(1), ...,a-i(j), fc+T, (a-i(j + 1), ...,a-i(fc))m(fcT2, ...,h)] 



We get: 



where (3j := 



11 



In the first term in this expansion of T2, we recognize — < a > P[i, ...,n — k], so that these terms sum 
up to — < i?0 > P[l, ...,n — k] when a runs over elements in Bk without regressions. In the second 
term (the sum over j < i), the same reasoning as for Ti shows that each term of the expansion is of the 
form < /3j;n — fc — 1 >, where Pj has no regression. Again, when a runs over elements in Bk without 
regressions, this provides (fc + 1)! such elements in the expansion. These terms are pairwise distinct and 
pairwise distinct from the elements showing up in the expansion of Ti, from which we conclude: 

T1-T2 =< i?§ > P[i,...,?r^]+ < i?^+^;n-fc + l > . 

The Theorem follows. 



6 A Magnus expansion for the evolution operator 

In the classical case, that is when the solution X{t) of a first order linear differential equation is obtained 
from its Picard series expansion, the resulting approximating series converges relatively slowly to the 
solution. This problem -let us call it the Magnus problem- is solved by reorganizing the series expansion, 
often by looking for an exponential expansion X{t) = expri(t) of the solution, known as its Magnus 
expansion. Many numerical techniques have been developed along this idea that go much beyond the 
formal-algebraic problem of deriving a formal expression for Q{t). However, deriving such an expression 
is a decisive step towards the understanding of the behavior of fl(t). This problem was solved, in the 
classical case, by Bialynicki-Birula, Mielnik and Plebanski O [T^ who obtained a formula for il{t) in 
terms of Solomon's elements D^. 

The purpose of the present section is to solve the Magnus problem for the analysis of solutions in 
time-dependent perturbation theory. This provides the general term of time-dependent coupled-cluster 
theory [26 . Our previous results pave the way toward the solution of the problem. Namely, as it appears 
from Thm l5.ll the natural object to look at is not so much the effective Hamiltonian 

n ^ lim PoHUGLie) 

or the effective adiabatic evolution operator Ugl, than the Picard- type series 

Notice that we define Pic as the sum of the < > over all the integers (and not over N*) in order to 
have the identity operator 1 — i/^o as the first term of the series. Of course, we have: 

UcLie) = P+{l-P){J2 < R% >)P = P+{1- P){J2 <R^>)P = P+ii- P) P^c P 

In other terms, we are interested in the expansion: 

UaUe) = P + (l - P)eM^e)P, 



where 



= log(^ < >) =< log(^ R^) > 



neN 

Since R;' * ... * R'^" = R^^Z^'-'lT+.-.+n.^i}^ ^ ^^^^^ expression of = log ^ follows: 



neN' SC[n-l] ' ' 

where one can recognize the hyperoctahedral analogue of Solomon's Eulerian idempotent [251 Chap. 3, 
Lem.3.14]: 

(_1)|5| 



SC[n-l] l*^' ' ' ^ 
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The analogy is not merely formal and follows from the isomorphism of Thm 14 . 1 1 together with the existence 
of a logarithmic expansion of soln, which is actually best understood from an Hopf algebraic point of 
view, see HOI [23 US [22: 

5^ sor = log(^ D^i). 

neN* riGN 

As a corollary of Thm 14. ll we also get the expansion of flu in the canonical basis of Q[i3„]: 



Proposition 6.1. We have. 



s 



riGW SC[n-l] 

(-1)1^1 



EE E 

neti' 5C[n-l] seB„,Reg{<y) = S 



\S\ 



The Proposition follows from the analogous expansion for soln |25j , together with the algebra isomor- 
phism Thmim 



soln = E E E 

new SC[ri-l] aeS„,Desc(cy)=S 



\S\ 



Corollary 6.2. The hyperoctahedral Magnus expansion of the effective Hamiltonian Ti. reads, when trun- 
cated at the third order: 

n = Ihn PHi{P+{l-P) exp(i7(i) +H^^ + ^[H^i2) +i^(12) +^(21) +^^(21) -^^(12) -i^(i2) -i^(2i) -if(2i)] + 

2 [H{123) + H(^i23) + i?(i32) + -ff(132) + -^(213) + -^"(213) + ^(231) + ^(231) + -^^(321) + -^"(321) + ^(312) + ^(312) 
+-ff(321) + -f^(321) + -ff(231) + -^"(231) + ^(123) + ^(123) + -^^(132) + -^"(132) + -^"(213) + -^"(213) + ^(312) + -^^(312)] 

[^(132)+-f^(132)+-^^(231)+-f^(231)+^(213)+-f^(213)+-f^(312)+^^(312)+-f^(132 
+^(2l3) + ^(312) + H(^l2) + i?(231) + -^(231) + -^"(321) + -^^(321) + ^^(123) + ^(123) + ^(321) + ^^(321)] ■ 
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